Relation of isospin-symmetry breaking correction for 
superallowed beta decay to energy of charge-exchange giant 

monopole resonance 

Vadim Rodin 1 

1 Institut fur Theoretische Physik der Universitdt Tubingen, 
Auf der Morgenstelle 14, D-72076 Tubingen, Germany 
(Dated: March 6, 2013) 

Abstract 

Having applied an analytical transformation, a new exact representation for the nuclear isospin- 
symmetry breaking correction 5c to superallowed beta decay is obtained. The correction is shown 
to be essentially reciprocal to the square of an energy parameter Qm which characterizes charge- 
exchange monopole strength distribution. The proportionality coefficient in this relation is de- 
termined by basic properties of the ground state of the even-even mother nucleus, and should be 
reliably calculated in any realistic nuclear model. Therefore, the single parameter Om contains all 
the information about the properties of excited + states needed to describe 5c- This parameter 
can possibly be determined experimentally by charge-exchange reactions. Basic quantities of inter- 
est are calculated within the isospin-consistent continuum random phase approximation, and the 
values of 5c are compared with the corresponding results by other approaches. 

PACS numbers: 24.80,+y, 23.40.Bw, 23.40.Hc, 21.60.-n, 21.10.Sf, 12.15.Hh 
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I. INTRODUCTION 



Superallowed + — > + (3 decays (SA (3 decays) allow to test fundamental properties of 
weak interaction, such as the Conserved Vector Current (CVC) hypothesis and the uni- 
tarity of the Cabibbo-Kobayashi-Maskawa (CKM) matrix (see, e.g., a recent review [l| 
by Towner and Hardy). The decay rates for SA Fermi transitions between T = 1 nuclear 
multiplet states have accurately been measured in a dozen of nuclei. Since the CVC hy- 
pothesis is only true in the isospin-symmetry limit, an uncertainty enters in the analysis 
of the experimental ft- values depending on a model calculation of the effect of the isospin 
breaking in nuclei. Although the breaking is weak, the current situation is such that the 
theoretical uncertainties in calculated correction terms predominate over the experimental 
uncertainties in the SA /3 decay data. This calls for a better accuracy of the theory applied 
to interpret the experimental results. 

Q ■ . 

From the 2009 survey of experimental data, Hardy and Towner [2J determined \V u d\ = 
0.97425 ± 0.00022, that, combined with the complimentary experimental data on \V US \ and 
\V U b\, gave for the norm of the first raw of the CKM matrix: 

\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 0.99990 ± 0.00060. (1) 

Thus, this test confirmed the unitarity of the CKM matrix with an accuracy of 0.06%. 

Isospin symmetry is slightly broken in nuclei, mainly by the Coulomb interaction. This 
leads to a small reduction of the nuclear matrix element Mp for SA Fermi transitions between 
the ground state (g.s.) of the even-even mother nucleus and its isobaric analog state (IAS) 
in the odd-odd daughter nucleus: 



\M F \ 2 



(2) 



where |M^,| 2 = 2T = \N — Z\ is the exact-symmetry value, with T = \T Z \ being the isospin 
of the g.s. of the even-even mother nucleus, and 8c > is the isospin-symmetry-breaking 
correction. 

There have been a number of methods used recently to calculate the correction 8c' shell 
model with Saxon- Woods and Hartree-Fock radial functions , relativistic Hartree and 
Hartree-Fock approach with the random phase approximation (RPA) [4], isovector monopole 
resonance (IVMR) model jlj], and self-consistent isospin- and angular-momentum-projected 
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nuclear density functional theory jg]. Still, there is a significant spread in the obtained values 
of 5c- Therefore, a better understanding of the aspects of nuclear structure important for 
more accurate evaluation of 5c is needed. 

The main purpose of this work is to derive a new exact representation for the correction 5c, 
which emphasizes the role of the physical charge-exchange monopole strength distributions, 
that can be probed experimentally. Having applied an exact analytical transformation, 5c 
is shown to be essentially reciprocal to the square of an energy parameter Qm which charac- 
terizes charge-exchange monopole strength distributions. The proportionality coefficient in 
this relation is determined by basic properties of the ground state of the even-even mother 
nucleus, and should be reliably calculated in any realistic nuclear model. Therefore, the 
single parameter Qm contains all the information about the properties of excited + states 
needed to describe 5c- The possibility of experimental determination of this parameter in 
charge-exchange reactions is discussed. Also in this paper basic quantities of interest are 
calculated for few nuclei within the isospin-consistent continuum-RPA, and the obtained 
values of 5c are compared with the corresponding results by other approaches. 



II. A NEW EXPRESSION FOR THE NUCLEAR COULOMB CORRECTION 5 C 

In all the experimental cases of interest, the IAS of the g.s. |0) of a even-even mother 
nucleus with T = 1 is an isolated low- lying state (or, in most of the cases, even the g.s.) in 
the daughter odd-odd nucleus. This fact allows for the following identical transformation of 
the representation for \Mf\ ■ 

\M F \ 2 = (0\f^\IAS)(IAS\f^\0) = 2T + S [ F +) - S { f\ (3) 
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Here, S ( p> = £ (j\T^\0) , Sf> = £' (i\T^\0) ,£' = £- £ runs over all physical 

j i i i i=IAS 

+ states, but the IAS, in the daughter nucleus, T~ = £ a r~ an d = (t^ ^ are the 
standard isospin lowering and raising operators, respectively. Note that Eq. ([3]) represents 
simply a version of the Ikeda sum rule for Fermi transitions: \Mp | 2 + — = 2T = 
N — Z. Hereafter we consider for definiteness the case N > Z, which can readily be 
generalized for the case N < Z. 

Therefore, the Coulomb correction 5c can be represented in the following form: 

^ = ^r(4~ ) -4 +) ) (4) 



It proves very useful to further transform Eq. as to explicitly relate 5c to isospin- 
breaking terms of the total nuclear Hamiltonian H, which include the Coulomb interaction 
and small isospin-violating part of the nuclear forces. For this, one introduces auxiliary 



operators V c 



(T) 
C 



± 



which are determined by these isospin-breaking terms, and 



uses an exact relation between the matrix elements of Vq and : 



(5) 



with u s = E s — E being the excitation energy of a state of the isobaric odd-odd daughter 
nucleus measured from the g.s. of the parent nucleus. The degree of fulfillment of Eq. (J5]) 
in a nuclear model can serve as an important check of the isospin-consistency of the model. 
In particular, Eq. (JSJ) ensures the equalities 
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where S^i] = £' 



(i\V c } |0) uf and S^l = £ 0'|Vc +) |0) are the energy- weighted 

j 

Coulomb sum rules. Therefore, one gets a representation alternative to Eq.(j3]): 
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2Tr 



C[-2] 



c(+) ^ 
D C[-2]s 



(7) 



Although physically Eq.(jl]) and Eq.([7]) are equivalent, it is preferable to use Eq.(J7]) in a 
model calculation as this representation is much less sensitive to possible residual isospin- 
inconsistences of the model. 

The charge-dependent isospin-breaking interaction is dominated by the Coulomb inter- 
action between protons. Because the Coulomb force is of the long range, the one-body 
Coulomb mean field is mainly determining the transition operators V c in Eq.((7]), which 
in this case also become one-body monopole charge-exchange operators. Therefore, most of 
the strength for the transition operators is exhausted by the corresponding giant isovector 
charge-exchange monopole resonance (IVMR), associated with the 2hu particle- hole exci- 
tations of proton-neutron type with J n = + . The importance of the IVMR as a doorway 
state for the isospin mixing of the IAS was realized already in the early years of the IAS 
studies jl, 8], and was re-emphasized recently in Ref. j^]. 

Further, one can introduce an auxiliary energy VIm defined as 
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It characterizes the charge-exchange monopole strength distributions in odd-odd iso- 
baric nuclei. Then the original expression for 5c © can be identically rewritten as 
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u\\M F \ 



and (IAS\V c ' } \0) = u A M F . From this expression one obtains \Mp 



Here, we again have used J2' = J2~ 

i i i=IAS 
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|0))), and, finally, arrives at the following expression: 



5, 



c 



O 2 
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\0)-u 2 A 



(9) 



Thus, one sees that 5c © is determined by two energies, Qm and ua, which are the only 
input in the problem related to the spectrum of + states in the odd-odd daughter nuclei, 
and by the properties of the g.s. of the mother nucleus via the g.s. expectation value of 



the commutator (0| V c + , V c |0). Both oja and Q,m can be determined experimentally 
(the former is already very accurately known; the value of the latter can be determined by 
charge-exchange reactions on the mother nucleus, see below). The numerical simulations 
(see below) indicate that a strong inequality Qm is fulfilled, that is to expect because 

of the high IVMR energy. 



Now we would like to evaluate the expectation value 



|0) in the dominant 



mean field approximation, Vc = J2 a Uc(r a )(l — r az )/2, with Uc{r) being the Coulomb mean 

field. The realistic potential Uc{r) resembles very much that of the uniformly charged 

Ze 2 

sphere, which is a quadratic function: Uc(r) = 7777(3 — (r/R) 2 ) inside a nucleus r < R, 

2,R 

where R is the nuclear radius. It turns out that if one extends this quadratic dependence 
also to the outer region r > R (instead of proportionality to 1/r), this gives numerically just 
a small deviation in the Coulomb sum rules. Thus, the Coulomb sum rules are determined 



by single-particle charge-exchange operators V c 



(=F) 

c 



U, 



C 



J2aUc( r a)T~a , where the term 
in Uc{r), proportional to r 2 , gives the dominant contribution to the sum rules. Further, one 
has 



C\ r a.Paz) 



and finally gets 



2T f 



<0| 



10) 



1 

2Tn 



U c (r)(Q n (r)-Q p (r))d 3 r = U c , 



(10) 



(11) 



where g n (r) and g p (r) are neutron and proton number densities, respectively, and the bar 
means averaging over the neutron excess density. 



The nominator in Eq. (Q is subject to a strong cancellation between the two terms. This 
can be best seen if one introduces the value of the Coulomb mean field averaged over the 
neutron excess density: 

By adding and subtracting Uc in the nominator of Eq. Q, one gets: 

Sc = 7^— 2 ((Uc - U~ c ) 2 + (U? - uj\)) . (13) 
Each of the two terms in the nominator of Eq. ( TT3l) is now much smaller than its counterpart 



in Eq. (£)]): (Uc — Uc) 2 <C U c as a consequence of the smoothness of the Coulomb mean 
field, and (Uc — u 2 A ) ^ °°a since Uc provides a leading contribution to the Coulomb 
displacement energy (see, e.g., {?], lsj|). The numerical simulations below show, that only 



for light nuclei the term (Uc — Uc) 2 is comparable with Uc — w A , for heavier systems the 



inequality U c — oj 2 a 3> \(U C — U c ) 2 \ holds. 

Different models must give similar results for U c and Uc provided that the basic nuclear 
geometry (like the mean radii of proton and neutron density distributions) can reasonably be 
reproduced. Therefore, different values of 5c obtained by different methods should mainly 
stem from differences in corresponding values of O^f. 

Returning to the question of the possible experimental determination of Qm, the charge- 
exchange IVMR was experimentally studied in ( 3 He,t) {J, ( 3 He,tp) 9- 11| and (t, 3 He) 11 ] 
charge-exchange reactions. The spin-flip IVMR was mainly excited in the experiments. 
Though the measurements are rather difficult, one may expect that the excitation of the 
non-spin-flip charge-exchange IVMR might be separated from its spin-flip counterpart (by 
means of polarized beams, or by comparing measurements at different projectile energies). 
Also, we note that the effective one-body transition operator leading to the IVMR excitation 
in the charge-exchange forward-scattering reactions is determined by the r 2 dependence of 
the Bessel function jo(r) 12], in accord with the r 2 dependence of Uc(r). In such a case, 
Qm can directly be obtained from the experimental cross-sections by the formula analogous 
to Eq. ©. 
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III. CALCULATION RESULTS 



We consider here by way of example four decays: 10 C— ?- 10 B, 38 K— )- 38 Ar, 66 As— > 66 Ge, and 
70 Br-)- 70 Se (two latter cases, both with A m 70, are taken to check the abrupt drop in 5c 
while going from 70 Br to 66 As as appeared in calculations of Ref. Q]). We use here semi- 
phenomenological nuclear mean field and apply the continuum-RPA with Landau-Migdal 
zero-range forces to calculate the quantities of interest: uja, U c ffTTl) . Uc ( 1T2|) . Sp, Sq^I, 
Sqm, Qm ©, and finally 5c- 

The mean field is chosen as described in Ref. [13J, and includes the fully phenomenolog- 



ical isoscalar part, with its parametrization tracing back to the Chepurnov's potential 14] . 
and both the symmetry potential and the mean Coulomb field calculated in the Hartree 
approximation. The chosen dimensionless intensity /' = 1.0 of the isovector part of the 
Landau-Migdal forces determines the symmetry potential via the isospin selfconsistency 



condition 



15]. Thus, the mean Coulomb field is the only source of the isospin breaking in 



the present model. 

The relevant equations of the genuine continuum-RPA are described in details in Ref. [13 



15] . Note, that we do not need any discretization of the s.p. continuum as done in Ref. 



because the equations are written in terms of the s.p. Green's functions. Since the nuclei 
in question are open-shell ones, one would in principle need to take into account the pairing 
correlations, and better use the continuum- QRPA [l6 18] instead of the continuum-RPA. 
However, the continuum- QRPA calculations are much more time consuming, and, more 
importantly, one can easily argue that the effect of nucleon pairing on the quantities in 
question must be small (since the pairing gap is much smaller than VLm)- Also, a moderner 



18], that allows for their A-dependence, 



choice of the isoscalar mean field parameters of Ref. 
is not expected to markedly affect the results. 

In Table [J the calculated excitation energies ua and charge radii r c are listed along with 
the corresponding experimental data. Table 1 also contains calculated values of U c (fTTTl and 
Uc (Il2i) . columns 6 and 7, respectively. The underestimate of the experimental IAS energy 



in the calculations for heavier nuclei reflects the Nolen-Schiffer anomaly 



19 



20]. 



The other calculated quantities of interest, which characterize the IVMR strength distri- 
butions, Sp , Sqmi, Sq^I, Qm ©, are listed in Table HT1 One can see a fairly good agreement 
between the corresponding entries in columns 2 and 4, and 3 and 5, respectively, in agree- 
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TABLE I: In columns 2 and 3 the experimental and calculated excitation energies uja, measured 
from the g.s. energy of the corresponding even-even nuclei in SA /3-decays from column 1, are listed 
(the calculated RPA values are corrected for the neutron-proton mass difference). The calculated 
charge radii r c of the even-even nuclei are given in column 5, and the only available experimental 
value for 38 Ar is given in column 4. In columns 6 and 7 the values of U c (jllj) and Uc (112p are 
listed, respectively. 

coa, MeV r c , fin L/g Uc~ 
exp calc exp calc MeV 2 MeV 



io c io B _j_ 397 _ 1 66 2.69 -8.96 -2.91 
38 K _^ 38 Ar 5 533 5 57 3 40 3 3Q 55 83 7 36 

66 As 66 Ge 9.609 8.97 3.93 142.80 11.91 
70 Br _^ 70 ge 10 .109 9.49 3.99 159.18 12.59 

ment with Eq. fl6]). This is a clear evidence of the isospin self-consistency of the applied 
continuum- RPA. Finally, the isospin-symmetry-breaking correction 5c is listed in columns 
9 (obtained directly from the RPA solution) and 10 (calculated from Eq. (TT3l) ). Both ways 
of calculating 5c agree well, again as a consequence of the isospin self-consistency of the 
applied continuum-RPA. 

TABLE II: Calculated Sp and S^_ 2 ^, in percents, are listed in columns 2,3 and 4,5, respectively. 
In columns 6,7 values of S^™ are given, and in column 8 the calculated energy Qm §B) is listed. 
The isospin-symmetry-breaking correction 5c, in percents, is listed in columns 9 (obtained directly 
from the RPA solution) and 10 (calculated from Eq. (|13p k 

Decay S F ,% 5 C[ _ 2] , % S c[0 ], MeV 2 Q M , MeV 5 C ,% 

+ + + rpa (in 



io c _^ io B Q65 35 o_o74 0.39 1.01 1.47 12.13 0.147 0.142 

38 K 38 Ar 2.I6 1.33 2.18 1.33 26.95 9.35 45.38 0.434 0.436 

66 As 66 Ge 7.38 5.44 7.43 5.43 96.51 19.77 61.90 0.992 1.007 

70 Br 70 ge 7 83 5 92 7 89 5 go 109 _ 47 21.17 66 . 66 .992 0.993 

Apart from the decay 66 As— > 66 Ge, the present results for 5c are close to those of Ref. 4| 



S 



calculated within the RH+RPA, and also are systematically smaller than those of Ref. j^]. 
As for the decay 66 As— > 66 Ge, the corresponding value of 5c is pretty close to the one for 
70 Br— >- 70 Se, in contrast to the conclusion of Ref. 141, but in a qualitative accord with a small 

n 

relative change in 5c between these decays as observed in Ref. p|. 

One can explain the difference between the shell model and the RPA results in terms of 
the difference in Qm- The former approach uses the differences in radial single-particle wave 
functions of the neutron and proton with the same quantum numbers, i.e. employs pure 
mean field picture. In this picture the collectivity of the IVSM is missing, and the effective 
Qm must be less than Qm of the continuum- RPA, in which a collective IVMR is formed 
by the repulsive residual particle-hole interaction and thereby shifted to higher excitation 
energy (see also similar arguments in Ref. js|). 

Note, that an estimate of the effect of the isospin splitting of the IVMR goes beyond 
the framework of the RPA. The splitting effectively pushes the monopole strength to higher 
excitation energies, and it is to expect that Qm will further slightly increase. This would 
lead to a corresponding decrease of 5c, bringing them closer to the estimates of Ref. j^. 



IV. CONCLUSIONS 



In present work a new exact representation for the correction 5c is derived, which em- 
phasizes the role of the physical charge-exchange monopole strength distributions. Having 
applied an exact analytical transformation, 5c is shown to be essentially reciprocal to the 
square of an energy parameter Qm which characterizes charge-exchange monopole strength 
distributions. The proportionality coefficient in this relation is determined by basic proper- 
ties of the ground state of the even-even mother nucleus, and should be reliably calculated in 
any realistic nuclear model. The possibility of experimental determination of the parameter 
Qm in charge-exchange reactions is discussed. Also in this paper basic quantities of interest 
are calculated for few nuclei within the isospin-consistent continuum-RPA, and the obtained 
values of 5c are compared with the corresponding results by other approaches. 



[1] I. S. Towner and J. C. Hardy, Rep. Prog. Phys. 73, 046301 (2010). 
[2] J. C. Hardy and I. S. Towner, Phys. Rev. C 79, 055502 (2009). 



9 



I. S. Towner and J. C. Hardy, Phys. Rev. C 77, 025501 (2008). 

H. Liang, N. Van Giai and J. Meng, Phys. Rev. C 79, 064316 (2009). 

N. Auerbach, Phys. Rev. C 79, 035502 (2009). 

W. Satula, J. Dobaczewski, W. Nazarewicz and M. Rafalski, Phys. Rev. Lett. 106, 132502 
(2011). 

N. Auerbach, J. Hiifner, A. K. Kerman and C. M. Shakin, Rev. Mod. Phys. 44, 48 (1972). 

N. Auerbach, Phys. Reps. 98, 273 (1983). 

R. G. T. Zegers et al, Phys. Rev. Lett. 84, 3779 (2000). 

R. G. T. Zegers et al, Phys. Rev. Lett. 90, 202501 (2003). 

K. Miki et al, J. Phys. Conf. Ser. 312, 092040 (2011). 

T.N. Taddeucci et al, Nucl. Phys. A469, 125 (1987). 

M. L. Gorelik, S. Shlomo and M. H. Urin, Phys. Rev. C 62, 044301 (2000). 

V.A. Chepurnov, Soviet J. of Nucl. Phys. 6, 696 (1967). 

M. L. Gorelik and M. H. Urin, Phys. Rev. C 63, 064312 (2001). 



V.A. Rodin and M.H. Urin, Phys. At. Nuclei 66, 2128 (2003), jnucj-th/0201065 
V. Rodin and A. Faessler, Phys. Rev. C 77, 025502 (2008). 

S. Y. Igashov, V. Rodin, A. Faessler and M. H. Urin, Phys. Rev. C 83, 044301 (2011). 
J. A. Nolen and J. P. Schiffer, Ann. Rev. Nucl. Part. Sci. 19, 471 (1969). 
S. Shlomo, Rep. Prog. Phys. 41, 957 (1978). 



10 



